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Abstract: We study a model describing a four-dimensional holographic supercon- 
ductor whose properties depend non-trivially on a particular spatial direction, namely 
a Josephson junction. We analyze the parameter dependence of the condensate and 
compare it to the condensed matter expectations. 
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1. Introduction 

Condensed matter systems often display a strongly coupled behavior near a quantum 
critical phase transition. The AdS/CFT correspondence [1, 2, 3] offers a powerful 
tool to study strongly coupled field theories. For reviews with particular emphasis 
on applications to condensed matter physics, see [4, 5, 6, 7]. The observation that 
a scalar field can break an abelian gauge symmetry near an AdS black hole horizon 
[8, 9] has led to the recent construction of a model of holographic superconducting 
phase transition [10]. 

The simplest model of the holographic superconductor minimally couples an 
abelian vector field and a massive scalar field to the Einstein-Hilbert action with 
a negative cosmological constant. In the pioneering works [10, 11], all the fields 
were supposed to depend on the radial AdS coordinate only, thus providing the 
gravity dual of a homogeneous superconductor. More recently, it has been shown 
that allowing the fields to non-trivially depend on the spatial coordinates of the 
boundary also gives rise to interesting phenomena in the dual field theory [12, 13], 
namely it makes it possible to describe a 2 + 1-dimensional holographic Josephson 
junction. The aim of the present paper is to extend their results to the 3 + 1- 
dimensional case and to show that the condensed matter expectations also match in 
this case. This allows us to further extend the number of phenomena which can be 
described via holography. 

The rest of the paper is organized as follows. In section 2 we write the gravi- 
tational Lagrangian, derive the set of equations of motion and review earlier results 
on the homogeneous solution to set the stage for the following discussion. Then, 
we move to the description of our inhomogeneous setting. In section 3 we present 
our numerical results for the superconducting phase and we compare them to the 
condensed matter expectations. We conclude in section 4. 

Note added: immediately after our work was completed, we became aware of 
[14], which has some overlap with our results. 
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2. The gravity model 



The minimal gravitational action describing a holographic model of a s-wave super- 
conductor contains an abelian gauge field and a charged minimally coupled scalar 
field 0. Allowing for a scalar mass term we then consider the action 



S = J d*x^r g [ JL (R+^- \F, V F^ - |D0| 2 - m 2 



(2-1) 



where the value of the cosmological constant has been set to 12, k 2 = 811G5 sets the 
five-dimensional Newton constant, 

= d^A v - d„A„ (2.2) 
= V M - iqA^ (2.3) 

are the gauge field strength and the covariant derivative respectively, and is the 
metric covariant derivative. The gauge field has been rescaled so that its kinetic term 
does not involve any constant and q is the charge of the scalar field in these units. 
The Einstein, Maxwell and scalar equations of motion reads 
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L 2 

V*%, = qi (0* - D 



(2.4) 



(D^Df, - m 2 ) = 

where G M „ is the Einstein tensor and is the energy-momentum tensor given by 

Tfiv = F W F/ - ^F pa F^ + 2D { ^D V) 4? - g, v \D<j>\ 2 - g^m 2 |0| 2 (2.5) 

When the fields are constrained to only depend upon the radial coordinate, the 
model (2.1) describes a homogeneous holographic superconductor. The common 
ansatz in that case is 

ds 2 = -f(r)e 2 ^dt 2 + + ^rfx 2 

/(r) L 2 

A t — h(r) Ai = i = r,x,y,z 



where the scalar can be taken to be real without any loss of generality and the 
function v vanishes if the Newton constant k — 0. There exist two phases separated 
by a critical temperature T c . When T > T c the scalar field vanishes (normal phase) 
and the analytic solution is 

}{r} ~L 2 V J r 2 



(2.7) 

h{r) = n 
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Below T c the scalar field acquires a non-trivial profile and the corresponding solution 
describes the superconducting phase. Because we are left with a set of non-linear 
coupled differential equations, there is not a general analytic solution, but many 
numeric reliable methods are available. We will take (2.6) to be one of our boundary 
conditions. 

We are interested in finding an inhomogeneous solution of (2.4), i.e. a solution 
in which the gauge and the scalar fields acquires a non-trivial dependence on the 
spatial coordinates. To simplify our task, we consider the theory in the probe limit, 
namely k — > 0; in that case the Einstein equations decouple from the gauge-scalar 
sector and we can consistently consider the plane-symmetric ansatz for the metric 

ds 2 = -f(r)dt 2 + -f- + T —dyi 2 (2.8) 
f(r) L 2 

We impose the boundary condition f{rn) = which determines, as a function of the 
horizon radius r# of the black hole, the explicit solution 



and the Hawking temperature 

T H = (2.10) 

which is identified with the dual field theory temperature. 

Without loss of generality we take the matter fields to depend on the spatial 
coordinate x but not on y and z. Thus, our ansatz reads 

Ai = Ai(r,x) i = t,r,x,y,z (f> = ip(r, x)e ie(r ' x) (2.11) 

In general, the phase 9 of the scalar field cannot be set to zero, but it disappears 
from the equations of motion once we shift the gauge field as A4 — > mi = Ai — did. 
Thus, a consistent subset of equations is given by only considering a non-vanishing 
m t which satisfies 

m'l + -m' t + i- d 2 x m t - ^- m t = (2.12) 

r r zj J 

where a prime indicates a derivative with respect to the radial coordinate r. Fur- 
thermore, the scalar equation reduces to 

„ /3 f'\ . L 2 o2 / m 2 m 2 \ 

r + {- r + j)^+^f d ^+{-Y + f)^ = ° (2 - 13) 

The solution to the system (2.12)-(2.13) is uniquely determined once we impose some 
boundary condition on the fields. We require that at large values for x the system 
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reduces to the homogeneous holographic superconductor described by (2.6). In that 
case, regularity imposes the gauge field m t to vanish at the horizon, while the scalar 
field has to satisfy 4>(th) = (/' \th) /m 2 ) <f>' \th)- At finite x, the general asymptotic 
solution near the boundary is 



m t {r, x) 



if)(r, x) 



fi(x) — 
^ + \x 



p{x) 



+ 



[x 



(2.14) 



-,4-A 



+ 



where m? = A(A — 4). According to the AdS/CFT dictionary, A is interpreted 
as the scaling dimension of the dual field theory operator which is acquiring a non- 
vanishing vacuum expectation value; thus, ip^ is interpreted as the condensate value 
(0a) of that operator and ip^ as its source. A spontaneous symmetry breaking on 
the boundary theory requires that we choose = 0. The field theory chemical 
potential and charge density are provided by fi(x) and p(x) respectively. They are not 
independent parameters: once we set one of the two, the other quantity is determined 
by the solution to the system (2.12)-(2.13). Thus, we set our last boundary condition 
by choosing the following function for the chemical potential [12] 



p(x) 



/'c 



1 



2Tanh [±] 













- ^Tanh 




— Tanh 
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)] 


a 


a 





(2.15) 



The peculiar feature of (2.15) is the following. For \x\ > / it approaches the constant 
value yUoo, while it quickly decreases at x ~ / to a constant value set by e. Thus, 
our model will describe two superconductors placed at large x and separated by an 
insulating barrier of size 2/ 1 . Although we do not expect our results to strongly 
depend on it, the function (2.15) is a well suited chemical potential to model an 
interesting class of condensed matter systems. 

The normal phase of our inhomogeneous solution is defined by the solution in 
which the scalar field vanishes identically. In this case, the boundary conditions we 
imposed only allow for the homogeneous solution 



h(r) = fi ( 1 



H 



(2.16) 



where [i is a constant. 

The superconducting phase corresponds to the hairy black hole solution, in which 
the scalar field assumes a nontrivial profile. When the temperature is lower than 
a certain critical value, the AdS black hole solution discussed above produces an 
instability because the scalar effective mass violates the BF bound. Thus, at the 
critical value for the temperature a phase transition occurs. 

^or low enough temperatures, the insulating barrier will switch to a superconducting one. We 
will work in a range of temperatures where this does not happen. 
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In the rest of this paper, we use the following scaling symmetry of both the bulk 
Lagrangian and the equations of motion 

r — > ar, t, x l — > at, ax 1 , L — > aL, q — )■ qj 'a,m — >■ raja; 
r — > br, t, x 1 — > t/fe, / — >• 6 2 /, m t — >• 6m 4 ; 

to set L = r# = 1 in performing numerics. In our approximations, the value of q 
does not affect our equations. 

For concreteness, we choose m 2 = —3. In this way, the dimension of the boundary 
operator which condenses is set to A = 3 2 and the value of the condensate is read 
from the asymptotic r~ 3 falloff of the scalar field, when a nontrivial solution exists. 



3. The superconducting phase 

The solution at a finite value of the spatial coordinate x strongly relies on the bound- 
ary solution at large x. The solution to the homogeneous solution is found by ex- 
panding the fields and the equations (2.12) and (2.13) near the horizon r = r# and 
by dropping the spatial dependence. Thus, we obtain 



m t {r) = m' t {r H ) 



(r -r H ) 2 (-6 + ?p(r H y 
th) H 



(r — ru 

I-—. — + 



3(r — r H ) | (r — r H ) 2 (33 — m' t {r H f 



4r H Q4r 2 H 



(3-1) 



+ 



We assign an initial value to m' t {rH) and VK r #); numerically integrate the equations 
from the horizon to the boundary and use a shooting method to seek a solution 
which satisfies = 0. Finally, we fit the asymptotic scalar solution to read off the 
operator condensate. The result is shown in figure 1. 

We solve the partial differential equations (2.12)-(2.13) by compactifying the 
radial and spatial coordinate as r' = 1 — 1/ r and x' = tanh [xj (4<r)) [12] and defining a 
grid on the resulting compact space; accordingly, we substitute derivatives with their 
discrete versions. At x — ±oo we impose that the fields approach the homogeneous 
solution discussed above, while at the boundary of the AdS space we assign to m t 
the values of the chemical potential. In the interior of the (r, x) space we use the 
relaxation method. Once we have the value of the solution at each point of the grid, 
we interpolate them and fit the asymptotic functions (2.14). 

An example of our solution is shown in figure 2. We draw a black line corre- 
sponding to the boundary chemical potential and plot the solution between x = — 10 
and x = 10, where it is indistinguishable from the homogeneous solution. In fact, it 
can be seen that immediately x is greater than /, the solution becomes very flat along 
that direction. That is a consequence of the shape of the chemical potential we have 
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Figure 1: The scalar condensate in the homogeneous solution as a function of the dimen- 
sionless temperature. 




Figure 2: An example of the gauge (left) and scalar (right) field solution. The black line 
in the left plot corresponds to the chemical potential. 



chosen as our boundary condition. A completely analogous argument also holds in 
the case of the scalar field. However, its shape distinguishes whether a non-trivial 
condensate appears in the dual field theory. For a sufficiently large I and low e in 
(2.15) its value at x = stays very close to zero 3 , resulting in a very low value 
for the condensate in that region. We show this phenomenon in figure 3. At x c± I 
the condensate starts to dip below the homogeneous limit and it quickly decreases 
reaching its minimum value at x — 0. Figure 3 shows two different regimes, namely 
/ = 1 and 1 = 6. In the former case the scalar condensate at x = is higher. Of 
course, we expected such a behavior. 

We can go a bit further and find a function which nicely interpolates our nu- 

2 we could also choose A = 1 in this case, but we do not expect our qualitative results to change 
with the dimension of the operator. 

3 the scalar field solution shown in figure 2 can be well approximated as <p(r, x) « 
4crArctanh 4 (x) ip(r), where (p(r) is the homogeneous solution. We did not find any similar in- 
terpolating function either for the gauge field or for the scalar field with I < 3. 
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Figure 3: The scalar condensate as a function of the spatial coordinate x for e = 0.6 and 
I = 1 (left) and I = 6 (right). 
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Figure 4: The scalar condensate as a function of the spatial coordinate x for I = 6 and 
e = 0.6. 



merical results at all Is. In the condensed matter literature, it is expected that the 
scalar condensate at x = and the length I are related by the simple relation 



(Os) 



T3 



Ae -i/m) 



(3.2) 



x=0 



where £ is the coherence length. Thus, we varied I from 1 to 6 using an interval 
of Al = 1/10 to test the relation (3.2) over a wide range of parameters. Figure 
4 shows the results of our scan. The points represent our numerical results, while 
the interpolating line represents the function (3.2) where the parameters assume the 
values A ~ 38.6 and £ ~ 0.88. Thus, we found a good agreement with the expected 
results. 



4. Conclusions 

We constructed a model of a four-dimensional inhomogeneous holographic supercon- 
ductor in the probe limit. Our Lagrangian was the same as used in the gravity dual 
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of the homogeneous holographic superconductor, showing that the same model can 
holographically describe a large class of condensed matter systems, both in 2 + 1 
and in 3 + 1 dimensions, by simply choosing appropriate boundary conditions for 
the fields in the ansatz. In the case of the boundary conditions we imposed, any 
anisotropics disappear at large values of the spatial coordinate the fields depend on. 
The lower dimensional counterpart was already studied in [12], and although our 
results are not directly comparable to theirs, we found many qualitative similarities. 

The critical temperature and the scalar condensate are not affected by the 
anisotropies we considered, because we took the solution to approach the homo- 
geneous one at large distances. Nevertheless, we have shown that an appropriate 
choice of the boundary conditions contains many parameters corresponding to the 
particular class of dual condensed matter systems we would like to model with our 
gravitational action. The full solution strongly depends on those parameters, allow- 
ing for the finding of some relation among them and the dual field theory data that 
the gauge/gravity duality predicts. Those relations have been shown to qualitatively 
agree with those expected from the condensed matter literature. This is, of course, 
a striking achievement. 

We did not attempt to study the conductivity and the way inhomogeneities man- 
ifest themselves there. As long as the the electric field is treated as a small external 
perturbation, it does not affect the superconducting background we constructed, but 
we expect it to depend on the spatial coordinates, even in the probe limit, and reduce 
to the conductivity of the homogeneous solution at large x. In particular, it would 
be interesting to look for the fate of the frequency gap over critical temperature ratio 
uj g /T c [15] in the interior of the x space and to give the electric gauge field a non- 
trivial dependence on another transverse spatial coordinate z and the wavevector k z 
to study the second order expansion of the Green function. Such coefficient contains 
the information relative to the magnetic permeability, and it can be shown that in 
the model we used here it contributes to achieving a negative refractive index. 

We only considered the simplest model in which a quadratic potential is given 
to the scalar field which represents the gravity dual of the boundary operator which 
condenses. Moreover, we only considered an abelian gauge field and the probe limit. 
The latter corresponds to give the bulk scalar field a large charge under the abelian 
gauge field. Both because going beyond these assumptions can display new inter- 
esting phenomena and because it seems difficult to derive a string theoretical model 
which reduces to this simple model, we think that it would be interesting to extend 
the present analysis to include more ingredients. 
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